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Abstract 

We apply an analytical approach for determining the near-field radiative heat transfer 
between a metallic nanosphere and a planar semi-infinite medium with some given surface 
structure. This approach is based on a perturbative expansion, and evaluated to first 
order in the surface profile. With the help of numerical results obtained for some simple 
model geometries we discuss typical signatures that should be obtainable with a near- 
field scanning thermal microscope operated in either constant-height or constant-distance 
mode. 
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I. INTRODUCTION 



Recent progress in the fabrication of a near-field scanning thermal microscope 
(NSThM) enables one to measure the radiative heat transfer between a cooled sample 
and a hot probe directly in the near-field regime, i.e. for distances in the nanometer 
rauge QQ. Due to thermally excited evaneseent waves. th.s regime one expects an 
energy transfer several orders of magnitude greater than the black body value [3|. For 
estimating the heat current in such a device several theoretical models are available, 
which describe the probe as a dielectric sphere, and the sample as a semi-infinite 



10| . These models can now be tested against 



dielectric body with a flat surface 
the data provided by the NSThM. 

In the literature, the near-field radiative heat transfer between a sphere and 
a structured surface has not been studied so far, although such a geometry is of 
considerable practical relevance and theoretical interest. Therefore, in the present 
paper we analyze the near-field radiative heat transfer between a spherical probe 
and planar samples with surface structures such as depicted in Fig. [H In particular, 
we discuss numerical results obtained for a planar surface structured by an infinite 
bar and a square pad, respectively. 

In this work, we use the general formulation of the near-field radiative heat trans- 
fer between a probe described as a spherical metallic nanoparticle within the dipole 
approximation and a second material as developed in Refs. js-iol, which is based on 
Rytov's fiuctuational electrodynamics JjJ. This formulation allows one to take the 
material's properties of the probe into account in terms of its electric and magnetic 
dipole moments, and the properties of the sample material, which is assumed to 
be a semi-infinite body with a given surface structure, through the local density of 
states 12| above that medium. Since the dipole moments of the probe or nanopar- 



ticle are known it remains to calculate the local density of states (LDOS) above 
the sample material. This is done within a perturbative approximation employing 
the Ewald-Oseen extinction theorem, as described in detail in Ref. [l^. It should 
be mentioned that our work is closely related to Ref. jl^, where the changes of 
linewidth and the lineshift for a molecule near a structured surface were calculated, 
since the electric LDOS can in principle be read off from the Green's function cal- 
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culated therein. However, in our work not only the electric but also the magnetic 
fields are determined, since it is known 
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15| that for metallic nanoparticles these 



magnetic fields can cause a heat transfer much greater than that due to the electric 
fields, as a result of the induction of Foucault's currents. 




a 



FIG. 1: Sketch of a probe-sample configuration with a square pad on a flat surface. 



This paper is organized as follows: In Sec. II we introduce the dipole model for 
the radiative heat transfer between a metallic nanoparticle and an arbitrary second 
material, and outline the strategy of the following calculations. In order to determine 
the electromagnetic LDOS, we deduce the appropriate dyadic Green's functions 
in Sec. Ill by means of a perturbative expansion due to Agarwal 13|, which we 
terminate after the first order. Within this approach, the LDOS itself is calculated 
and numerically evaluated for different surface profiles in Sec. IV, where we also 
deduce tentative criteria justifying the restriction to the lowest-order contributions. 
In Sec. V we then present some numerical results for the near-field heat transfer 
between a metallic nanosphere and a structured surface, and predict signatures that 
should be observable with a NSThM operated in either const ant- height or constant- 
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distance mode. 



II. HEAT TRANSFER BETWEEN A NANOPARTICLE AND A SECOND 
MEDIUM 



It has recently been shown 



Isijiol 



that the near-field radiative heat transfer be- 
tween a metallic nanoparticle with radius R and temperature Tp and a dielectric 
material with temperature Ts can be described within the dipole model for particle- 
sample distances a ^ R. In this model the electricpolarisability a-p and the mag- 



netic polarisability fip of the particle are given by 16 1 



ap(c.,fl) = 4rf^^lM^ d) 

epiu) + z 



and 



R' 

/ip(w,i?) = — - 



3^-^ + 3^^ cot 



R Kd'Juj 



s\ 



(2) 



with the permittivity of the probe ep{u) and d'^ico) := c{u ep{u) — 1)"^; as usual, 
c denotes the velocity of light in vacuum. Assuming that the two bodies are in local 
thermal equilibrium, the mean energy rate P fiowing from the hot to the cold body 
is given by the relation 

POO 

P= du[Q{u,Tp)-Q{uj,Ts)]2uj{a'p{u,R)D^{u,r) + iJ,'^{u,R)D^{u,r)). (3) 
Jo 

Here, the sign of the mean energy rate P is determined by the difference of the Bose- 
Einstein functions Q{uj, T) := hu{exp{hw f3) — 1)"^ with the temperature := ksT, 
so that the sign is positive for Tp > Tg and negative otherwise. The material's prop- 
erties of the nanoparticle are taken into account in Eq. ([3]) by means of the imaginary 
part (as indicated by the double prime) of the electric and magnetic polarisabilities, 
ap{u,R) and ^p{u,R), whereas the material's and geometrical properties of the 
second medium enter into this expression via the electric and magnetic local density 
of states D^{uj, r) and D^{uj, r) at the point r above this medium, where the probe 
is located. Therefore, the thermal radiative heat transfer between a nanoparticle 
and an arbitrary second medium can be calculated if the electromagnetic local den- 
sity of states above that medium is known. We point out that in the opposite limit 
a R the dipole model discussed here is not valid, whereas the so called "proximity 



approximation" should prove to be useful then [21 1. 



Now the local density of states can be calculated with the help of the relations 

G^(r,r')G'^(r,r')^ 
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and 



DE(^^r) = ^e'g'H /dVTr 
vrc^ J 

u e'g'(a;) 



d r Tr 



V X G^(r,r')(V x G^(r,r'))^ 



(4) 



(5) 



in the case of local equilibrium inside a heat-radiating local medium surrounded 
by vacuum. Here /io is the magnetic permeability of the vacuum, and 63(0;) is the 
relative permittivity of the material considered. These relations can in principle 
be evaluated by determining the electric dyadic Green's function G^(r, r') with the 
source points r' inside the medium and the observation point r outside the medium, 
implementing the tensor product with its hermitian conjugate and integrating over 
the volume of the medium. Since, here only a local equilibrium inside the medium 
is assumed, one can use these relations to investigate for example the heat transfer 
between bodies kept at different temperatures. 

Considering a medium surrounded by vacuum, the non-equilibrium expressions 
in Eqs. (jlj) and can be decomposed into an evanescent and a propagating part 



2l|. This evanescent 



by transforming the volume integral into a surface integral 
part of the local density of states coincides with the equilibrium expression, since 
the evanescent modes are bound to the surface of the medium and are therefore not 
relevant for preserving a global equilibrium situation. Thus, the evanescent part 
of the local density of states above a material surrounded by vacuum can also be 



calculated by means of the equilibrium expressions 
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D^{uj,r) = — ImTrG|^(r,r) 



and 



vrc- 



^ImTrGf(r,r) 



(6) 



(7) 



which state that the electric and magnetic local density of states at the point r 
is given by the imaginary part of the trace of the renormalised Green's functions 
G^(r, r) and Gj|(r,r), where the renormalisation procedure is defined as 201] 



G 



E/H 
R ' 



r, r 



lim 



G^/»( 



r, r 



G, 



E/H 



r, r 
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with GQ^^(r, r') denoting the Green's function of the vacuum. Since the Green's 
function with the observation and source point located above the medium con- 
sists of an incident and a reflected part, with the incident part coinciding with 

E /H 

Gg (i",!"')) the renormalised dyadic Green's function coincides with the reflected 
Green's function, so that the index "R" can be understood as both "renormalised" 
and "reflected". 

Obviously, the relations in Eqs. and ([7]) are much easier to evaluate than the 
expressions in Eqs. and (E]), giving reliable results for such distances above the 
material, at which the evanescent modes dominate the local density of states. Since 
we are especially interested in the evanescent regime, we focus on the equilibrium 
relations in order to determine the radiative heat transfer between a nanoparticle 
and a structured surface, keeping in mind that the results hold in the evanescent 
regime only. Thus, it is necessary to calculate the reflected electric and magnetic 
Green's function with observation points r and source points r' located outside the 
material of interest. For an electric source current jeo and a magnetic source current 
jmo located at the point r', the reflected fields and reflected dyadic Green's functions 
are connected by the relations 

Ek{co, r) = iw/ioG^(r, r')jeo (9) 

and 

Hr(w, r) = iufioGlir, r')imo- (10) 

Furthermore, in vacuum the electric and magnetic dyadic Green's functions are 



related by 



23| 



G«(r,r') = -^VxG^(r,r')x V (ll) 



with ko := u/c. 

In the following, we will therefore perturbatively evaluate the reflected electric 
field generated by an electric current jeo located at the source point r' above a semi- 
infinite medium with a structured surface, determine the electric dyadic Green's 
function by means of Eq. and from this result deduce the magnetic dyadic 
Green's function from Eq. (fTTl) . Finally, we calculate the local densities of states 
above the structured surface with Eqs. ([6]) and (IT!). These densities, in their turn. 
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then allow us to determine the near-field radiative heat transfer between a metallic 
nanoparticle and a semi-ininite medium with a structured surface. 



III. GREEN'S FUNCTION ABOVE A STRUCTURED SURFACE 

In order to calculate the electromagnetic fields above a structured surface, we 
assume that the surface profile is given by an expression hf{x,y), where the dimen- 
sionless function f{x,y) varies between zero and unity, and h is the characteristic 
scale of the profile variation, as sketched in Fig. O Moreover, we assume that h be 
small compared to all other relevant length scales of the problem, so that we can 
apply a perturbation expansion. This approach has been worked out in great detail 



by Agarwal 13|, so that it suffices here to mention only those elements that are 



indispensable to follow our line of reasoning. 
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Within this approach the so-called Ewald-Oseen extinction theorem p/7 
employed, allowing us to restate the boundary conditions of the electromagnetic 
fields in the given geometry as integral equations. For the case of a non-magnetic, 
isotropic, local, and linear material this theorem states that for observation points 
r outside that material, i.e. r ^ V (see Fig. [2]), one has 



E(r) = Ei(r) + ER(r) 



1 

^0 JdV 



Ei + -2 V X V X / dS' 

Kr 



ET(r) — ^(r-r ^ 



dn' dn' 



(12) 



Here, the field E(r) outside the medium is simply the sum of the incident field 
Ei(r) and the refiected field ER(r). The latter is described by the surface integral 
in Eq. ( fT2l) . so that the reflected field ER(r) can be calculated by means of the free 
Green's function 

ifco|r-r' 

g{r-r') = ^ (13) 

47r|r — r'l 

and the transmitted field ET(r). In addition, ^ symbolizes the normal derivative, 
taken in the direction of the unit normal of the surface, 

'i + A2|v;i/(x',t/')P' 



with Vji := {dx>,dyr,Oy and the unit vector in 2;-direction. 
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FIG. 2: Sketch of an example structure, showing a particular value 
of the surface profile hf{x, y) at x = xq and y = 0. 

On the other hand, for observation points r within the medium, i.e. r & V, the 
Ewald-Oseen theorem gives 



= Ei + -^V X V X 



d^' 



dV 



Exfr' 



dg{r 



.9ET(r') 



(15) 



dn' dn' 
In fact, the surface integrals in Eqs. f lT2|) and f|T5l) have the same structure, but one 
has to keep in mind that r ^ in Eq. 0121) . whereas r G ^ in Eq. fllSp . so that both 
integrals give different results. By means of Eq. ffT^ the transmitted field ET(r) can 
be computed if the incident field Ei(r) is given, so that the sought-after reflected 
field ER(r) can then be determined with the surface integral in Eq. ( |T2i) . 

Now the transmitted and reflected fields in Eqs. (fT2l) and ( |T5ll are expanded in a 
power series 



E 



T/Rl, 



(16) 



n=0 



Furthermore, all quantities in the Ewald-Oseen extinction theorem are expanded in 
a Taylor series with respect to the small quantity hf{x,y), giving 



ET(r') = Ei?\x',y',z' = 0) - hfix',y')^{x',y',z' = 0) + Oih' 
for the transmitted field and 



dE 



(0) 
T 



dz' 



(17) 



9{r 



-hf{x',y')^9{r-r') 

z'=0 OZ' 



for the Green's function. The normal derivative times the surface element dS' yields 



d^^ = -dx'dy' [e, + V[hf{x', y')] ■ V 



(19) 



Substituting these expressions into Eqs. ffT^ and ffT^ . and comparing coefficients 
yields constitutive equations for the fields in different orders of h. 
Next, it is useful to expand the fields in plane waves according to 



ET(r) = 
ER(r) = 
Ei(r) = 

and to utilize the Weyl expansion 
^(r - r') = 



(27r)2 



(20) 



d^K i 



^i(K,-{x-x')+fc2ok-z'|) 



(21) 



(27r)2 2ho 

for the Green's function, having introduced the notation k'^Q = k'^—H?, k^ = /Cge^— k^, 
X := {x,yY, and k, = {kr^kyY. It has to be emphasized that the use of these 
expansions down to the surface relies on the Rayleigh hypothesis, as discussed in 
Ref. 0. 

With the help of the plane-wave expansions the relations for the Fourier compo- 
nents of the refiected fields in terms of the incident field are easily calculated; for 



details we refer again to the work by Agarwal 
obtains the expression 



13|. For the zeroth-order field one 



kz — kzo 2kzoikz — kzo) , 

1 + ^ ( K - kzkz^^z] 



kz + kzo kQ^kzoSr + k 
and for the first-order field 



(22) 



4^(-)=i(^s-l)/|^F 



K, — K,') Lf/t) 



2k' 



zO 



-1 + 



2k' 



zO 



(23) 



with K' := (k', fc^)*, Kq := (k', /c^q)*; and the dyadic operator 
1 



kzQ^s ~ kz 



[k, + kzkzo){t — Bz ® e^) — K ® K + kzBz ® K + ErK Gz ® ^z 



ErkzoK ® Gz 



(24) 



Here ® symbolizes the dyadic product of two vectors, and 1 the unit dyad. Of course, 
also higher orders can be calculated within this approach [l3|, but the higher-order 
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contributions become increasingly cumbersome. Thus, this approach is particularly 
useful if meaningful results can already be obtained to first order. Hence, we re- 
strict ourselves here to these first-order fields, and try to give approximate criteria 
justifying this termination of the series later on. 

From Eq. ([9]) it is clear that the Fourier component of the electric dyadic Green's 
function for the reflected fields given in Eqs. fl22l) and fl23|l can be read off if we 
consider the incident electric fields S\ generated by a delta-like source current jeo 
located at r" and put the Fourier component of the result into Eqs. fl221) and fl2^ . 
According to Eq. ([9]), this field can be stated as 

Ei(a;,r) = ia;/ioG^(r,r")-jeo, (25) 

where here Gq is the free dyadic Green's function, which can directly be obtained 
from the relation 

g{Y-Y"). (26) 



G^(r,r")= (l 



IP- 



Inserting the Weyl expansion from Eq. (|2T1) for the Green's function g{r — r"), such 
that z" < z, yields 

G^(r,r") = y ^ (e^®e^ + e|,(-A:,o)®e||(-A:,o)). (27) 

Here we have defined the unit vectors in vacuum for the TE- and TM-modes as 

e^:=-{ky, -K,OY (28) 

and 

eii{k):=-^^{Kk,kyk,K^y. (29) 

Using this expression for the free Green's function in Eq. ( l25l) allows us to identify 
the Fourier component of the incident field, which is 

£i{k) = -ufio —, (e_L ® + ei\{-kzo) (g) ei\{-kzo)) ■ jeo- (30) 

//C20 ^ ^ 

Now, it is a straightforward exercise to calculate the reflected fields. Substituting 
Eq. ( I30l) into Eq. ( l22l) gives the zeroth-order field 

= -^/"o ( ^±e_L ® e_L ryen {k^^ ® ey [-k-^) ) ■ jeo, (31) 
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where we have introduced the usual Fresnel reflection coefficients for the TE- and 



TM-modes, defined as 



and 

^s^zo — kz 



" ^skzo + kz ^ ^ 

The ffist-order field can be calculated by substituting Eq. (|30!) in Eq. (123]) . giving 

J (27r)^ 2kzo (34) 

X (t'^el ® el + t\e\{-k',) ® e\\{-k',)) ■ jeo, 

with the transmission coefficients for the TE- and the TM-modes defined as 

kzO "I" 



t|| := , ■ (36) 



and 



Due to the definition of the vector ey for the TM-modes in Eq. ( 1291) the TM-mode 
transmission coefficient defined here does not coincide with the standard formulation 
of the transmission coefficient. Finally, the reflected electric dyadic Green's function 
can be read off, giving the zeroth-order expression 

G^'°(r,r") = (rxex®ej. + r,|e||(M®e|,(-A;,o)), (37) 

and the ffist-order expression 

}2,^ r a2,j 



g'^(r, r") = (1 - ^s) / (1;^ F(A^ - ^^0 H>^') -l-eH--'=-^)e-('^'-+'=-^) 

X (t'^e'^ ® e'^ + t\e\{-k'^) ® G\\{-k',)) . 

(38) 

Since the magnetic dyadic Green's function is linked with the electric one by means 
of Eq. ( ITTl) . we now have all ingredients to calculate the LDOS above a structured 
surface up to ffist order. 
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IV. LOCAL DENSITY OF STATES ABOVE A STRUCTURE SURFACE 



The zeroth- and first-order Green's functions in Eqs. fl37l) and flHHjl can now be 
used to calculate the electric and magnetic local density of states 

D^/«(u;,r)^Z}o''/^(u;,r) + Z}f/«(u;,r) 



Im Tr G|/^'° + h^lm Tr gI^^'^ 



(39) 



R 



by means of the equilibrium relations in Eqs. (E]) and (I7j). As expected, and as a 
confirmation of the validity of the approach, to zeroth order we obtain the well- 
known expressions [l^ 



and 



Att'^c'- 



■Im / dk. 



Im / dk. 



A;^e-2T" 
27 

fc^e-^T"^ 
27 



2/C /Cn 



rii -|- 



2fex ^0 



(40) 



(41) 



for the LDOS at the distance a = —z above a semi infinite body, with 7 = \/k'^ — k^. 
The first-order contributions are 



Df{u,a) =h\m 



1-eg fd'n fd^n' 



47r2 / 47r2 



F(h. - At') eK'^-'^'y^eKk.o+Ko)'^ 



kzkzo{i'i + k^k^o) + kJt^^i^K + k^k^o) 



(42) 



E/iIm 



l2, / 



1-es /"d^K [d^'K 



Fin - K')e'^''~'^'^-''hiK, k' 



and 



(w,a) =hlm 



^2, J 



es - 1 /"d^K fd^K' 



Fin — n') e^('^~'*')'^e''''^'=°"'"'^^o^° 



2^20 2^20 



47r2 7 47r2 



1^ 



^0 



+ k'^^ [k, - kzo) ^ + K [kz - kzo) 



./2 



(43) 



— esk^K ■ k! 



E/ilm 



es - 1 /"d^K /"d^K' 



47r2 / 47r2 



F{n- n')Q'^''-'^'^-''l^{K, k') 
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From these two equations (1421) and fl45|) it can be seen that in first order the electric 
and magnetic local density of states are given by an integral of the Fourier transform 

F(k-k') = Jdx Jdyf{x,y)e-''''<'^~'^'^ (44) 

of the profile f{x,y), multiplied by the exponential factor exp(i(K; — k') ■ x) and the 
function or respectively. Due to the factor exp(i(A;20 + Ko)^)^ 

these two functions are exponentially damped for k and k' with modulus much 
greater than the inverse observation distance in the near field, since in the 
evanescent regime kzo ~ in and k'^Q ~ in'. Therefore, the smaller the observation 
distance the more Fourier components contribute to the integral. Hence, from the 
structure of the first-order integrals it can be expected that for distances smaller than 
the characteristic width of the surface profile function the electric and magnetic local 
density of states resembles the surface profile. 

This behaviour is confirmed in Fig. [3], where the numerical evaluation of the 
electric density of states from Eq. ( l39l) above a square pad with a height h = 5 nm and 
width w = 15 nm on a plane surface is plotted. Here, we have used the frequency u = 
10^^ s~^ and the Drude permittivity es of gold. Obviously, at an observation distance 
of a = 40 nm the values of the local density of states resemble a two dimensional 
bell-shaped function, whereas at an observation distance of a = 5.5 nm, being much 
smaller than the width of the square pad, the values of the local density mimic the 
underlying structure, apart from softened edges. 

In the following, we discuss the magnitude of the contributions of the zeroth and 
first order of the electric and magnetic local density of states. For that purpose, we 
assume an infinitely extended bar with a width w; on a plane surface modelled by 
the profile function 

exp{d{\x\ - f )) + 1 

assuming d = 10^ m~^ and w = 30 nm. The height of the bar (see Fig. H]) is chosen 
to he h = 5 nm. 

The numerical values of the local density of states at a constant observation dis- 
tance a = 10 nm above the base plane are shown in Fig. [5l using again the frequency 
u = lO^^s"^ and the Drude permittivity for gold. Firstly, one observes that the 
values of the magnetic LDOS are much greater than the values of its electric coun- 



13 




(a) Observation distance a = 40 nm (b) Observation distance a = 5.5 nm 

FIG. 3: Electric LDOS from eq. (j39p above a square pad with an edge length of 15 nm and 
a height of /i = 5nm for two different observation distances a from the base plane. The 
frequency considered is w = lO^^s"^. 

profile 




-20 20 

lateral position x (nm) 



FIG. 4: Plot of hf{x,y) with the profile function defined by Eq. (j45|) with d = 10^ m ^, 
w; = 30 nm, and h = 5 nm. 

terpart, which is typical for metals, whereas for a polar dielectric bodies the electric 
LDOS usually gives greater values than the magnetic one. Secondly, at the given 
observation distance we find values for the first-order contribution to the electric 
LDOS, which are of the same order of magnitude as the zeroth-order contribution. 
On the other hand, the values of the first-order magnetic LDOS are significantly 
smaller than the zeroth-order values for all lateral positions. Thirdly, the values of 
Df and give an equally good image of the bar on the plane surface, the width 
of the two bell-shaped curves being approximately the same. Therefore, in this case 
both first-order contributions give qualitatively similar results. 

In Fig. |6] we plot the ratio of the first- and zeroth-order LDOS Di/Dq for the 
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-20 20 -20 20 

lateral position x (nm) lateral position x (nm) 

(a)Electrical LDOS (b)Magnetical LDOS 

FIG. 5: LDOS of an gold half space structured by a gold bar (see Fig. H]). The frequency 
used is Lo = 10^^ s~^, and the observation distance 10 nm. Observe the different scales for 
the electric LDOS (a) and the magnetic LDOS (b) 

same surface profile directly above the bar, i.e. at a; = 0, for observation distances 
ranging from 6 nm to 100 nm. As expected, this ratio increases for decreasing obser- 
vation distance, suggesting that for small distances higher-order terms have to be 
considered. These higher-order terms can in principle be calculated by an iterative 



scheme deduced by Greffet [19|. We suggest that it might be sufficient to consider 
only Do and Di as long as the ratio Di/Dq does not exceed the ten percent level. 
This means that the numerical results depicted in Figs. O a) and b) refer to distances 
where higher-order terms should be taken into account. Furthermore, it is evident 
from Fig. |6] that for all distances the ratio of the leading two contributions to the 
magnetic LDOS gives much smaller values than the corresponding ratio for the elec- 
tric LDOS. Therefore, it can be concluded that the approximation committed when 
considering only the zeroth- and first-order terms holds for the magnetic part for 
much smaller distances than for the electric part. This also means that the under- 
lying structure becomes important in the electric LDOS for much greater distances 
than in the magnetic LDOS. 

So far we have determined the LDOS at a constant height above the base plane of 
the structured surface. A very common operation mode in near-field microscopy is 
the constant-distance mode, where the separation between the tip and the individual 
features of the sample is kept constant 2[. In order to calculate the LDOS relevant 
for this mode, we use the observation distance a + hf{x,y) instead of a = const, in 
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0.01 I ' ' ' 1 

6 12 18 

a/h 

FIG. 6: Dependence of the ratio of the first and zeroth order of the LDOS on the 
observation distance a calculated above the profile shown in Fig. H] for x = and 
frequency = 10^^ s~^. 

Eq. Fig. [7] shows the LDOS obtained for this constant-distance mode above 

the structured half space depicted in Fig. HI 
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FIG. 7: Electric local density of states {D^: dash-dotted line; Df: dashed line; 
Dq + Df: solid line) above the structured surface as depicted in Fig. d] at a 
frequency ofuj = 10^^ s~^ and a constant separation of 10 nm between observation 
point and surface. 

In the constant-height mode (see Fig. [5] (a) and (b)) the first-order term Di gives 
a rough image of the underlying surface structure, while the zeroth-order term Dq 
is constant. However, in the constant-distance mode (see Fig. [7]) the qualitative 
behaviour of the two contributions is more complex, due to the variation of the 
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observation distance. Here, the zeroth-order term Dq coincides more or less with 
the inverse of the underlying structure, and the first-order term Di gives relatively 
large values near the edge of the bar, due to the variation in observation distance. 
Therefore, at least two regimes of observation distances can be distinguished: At 
large distances, where Di/Dq is small, the LDOS is dominated by the zeroth-order 
term, giving values which coincide approximately with the inverse of the surface 
structure. At distances where at least the first-order term has to be taken into 
account, the variation in distance leads to a rather complex pattern of the LDOS. 



V. CONSEQUENCES FOR THE NEAR-FIELD SCANNING THERMAL 
MICROSCOPE 



Up to this point, only the LDOS at a frequency corresponding to 300 K have been 
calculated. Before we discuss some numerical results for the full near-field radiative 
heat transfer between a metallic nano sphere and a structured surface, we emphasize 
the restrictions of the model: Firstly, due to the dipole approximation this model 
is only valid for distances a significantly greater than the radius of the sphere R. 
Secondly, for distances and sphere radii smaller than the mean free path of the 
conduction electrons, nonlocal and quantum mechanical effects become important 
and have to be implemented. Thirdly, due to the perturbative approach the results 
apply under the condition that the height h of the profile is the smallest length scale, 
i.e. h ^ min{a, Ath, w} with the width of the surface structure w and the thermal 
wavelength Ath ~ ^/3c. In particular this means that the pertubative approach is 
not valid for distances a <^ R; in this limit the "proximity approximation" can be 



utilised 



21| . Nonetheless, it is reasonable to explore the perturbative predictions in 



some detail. 

To this end, we evaluate the near-field radiative heat transfer between a gold 
nanoparticle and a gold sample as given by Eq. ([H]), using the polarisabilities of the 
sphere ([T]) and with the Drude permittivities ep and es for gold. As surface profile 
we employ the infinitely extended bar from Eq. ( H5l) with parameters as in Fig. HI 
In order to relate the results for the energy flow P to the numerical results for the 
LDOS we evaluate Eq. (jS]) for a constant height, a = 10 nm, and a constant distance, 
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a = 10 nm + hf{x,y), where a nanoparticle radius of i? = 10 nm has been chosen. 
As shown in Fig. [HI in both cases the near-field radiative heat transfer is dominated 
by the magnetic contribution, this being a typical feature for good metals like gold 
(at the given distance). Apart from this the curves displayed in Fig. |H]refiect the 
corresponding plots of the LDOS at a frequency near the thermal frequency Uth- This 
is a consequence of Eq. ([3]), since the main contributions to the frequency integral 
stem from frequencies near the thermal frequency as long as there are no resonances 
in the thermally accessible frequency regime. We also implemented the "proximity 
approximation" for the given geometry numerically (with the same parameters as 
in Fig. IH]), obtaining qualitatively similar results for the near- field radiative heat 
transfer as those given by the dipole model. 




-20 20 -20 20 

lateral position x (nm) lateral position x (nm) 

(a)Heat transfer in the constant-height mode. (b)Heat transfer in the constant-distance mode. 

FIG. 8: Near-field radiative heat transfer between the structured sample depicted in Fig. H] 
(Ts = 100 K) and a spherical probe (Tp = 300 K) with a radius of 10 nm. The height and 
the separation was 10 nm in each case. 

Hence, the conclusions deduced in the last section for the local density of states 
can be transfered to the near-field radiative heat transfer between a sphere and a 
structured surface, modelling a NSThM tip and a structured sample. Thus, for ex- 
ample we conclude that the lateral resolution of a NSThM increases for decreasing 
distances, as long as this microscope can be described within the dipole model. The 
lateral resolution is comparably good for the electric part Pe and the magnetic part 
Ph, with the topology of the surface profile being well resolvable for distances much 
smaller than the width of the surface pattern. On the other hand, for non-metallic 
materials, i.e. in a situation with Pg ^ Ph, the underlying structure becomes im- 
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portant for much greater distances than for metalhc materials, as can be concluded 
from Fig. [HI 

Furthermore, a NSThM gives very different images of the underlying structure 
when either the constant-height or the constant-distance mode is used. In partic- 
ular, for distances where Dq ^ Di the values for P should resemble the inverse 
surface structure in const ant- distance mode, whereas for distances where Di be- 
comes important, the measured signal can be rather complex in that mode. On 
the other hand, in constant-height mode the signal should approximately resemble 
the topology of the underlying structure at distances where Di becomes important. 
These are strong predictions that are amenable to immediate verification in current 

n 

experiments [2|, and may help to correctly interprete the various types of signals 



obtainable with a NSThM 
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